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Abstract 



With only three flavors it is possible to account for various neutrino os- 
cillation experiments. The masses and mixing angles for three neutrinos can 
be determined from the available experimental data on neutrino oscillation 
and from the astrophysical arguments. We have shown here that such masses 
and mixing angles which can explain atmospheric neutrino anomaly, LSND 
result and the solar neutrino experimental data, can be reconciled with the 
i?-parity violating Supersymmetric Models through lepton number violating 
interactions. We have estimated the order of magnitude for some lepton 
number violating couplings. Our analysis indicates that the lepton number 
violation is likely to be observed in near future experiments. From the data 
on neutrino oscillation and the electric dipole moment of electron, under some 
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circumstances it is possible to obtain constraint on the complex phase of some 
supersymmetry breaking parameters in i?-parity violating Supersymmetric 
models. 

PACS number(s): 14.60.St, 14.60.Pq, 12.60 Jv 
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I. INTRODUCTION 



Although in the Standard Model of electroweak and strong interactions the neutrinos 
are massless to all orders in perturbation theory, in its extension, the neutrinos may ac- 
quire small masses with see-saw type mechanism in presence of sterile neutrinos. Also such 
masses can be present in the minimal supersymmetric model with the renormalizable lepton 
number violating terms in the lagrangian. On the other hand, the astrophysical and cos- 
mological considerations also strongly suggest the existence of massive neutrinos. Presently, 
there are some possible evidences PJ of massive neutrinos and the mixing of different flavor 
of neutrinos particularly coming from the anomalies observed in the solar neutrino flux |§, 
in the atmospheric neutrino production H[| and in the neutrino beams from accelerators 
and reactors ||. Although some of the evidences like those coming from solar neutrinos 
and accelerator data has been explained || considering one massive and two nearly massless 
neutrinos but it is in general difficult to fit various neutrino data considering three neutrinos 
as particularly the first three evidences are best fitted by three different mass gaps for neu- 
trinos. However, the conventional approach to analyse various observed neutrino anomalies 
in the experiments is to parametrise those in terms of oscillation of two neutrino states only. 
This assumption may not hold good while fitting several observed anomalies simultaneously 
and consistent three flavor mixing scheme [7| for three neutrinos to analyse various data 
is essential. Several authors [j§HlOH have tried to fit various experimental data on neutrino 
oscillation in the three flavor mixing scheme. It is interesting to note that including the 
recent CHOOZ and the SuperKamiokande || result on neutrino oscillation alongwith 
other experiments in this direction, it is possible to find the mass square differences and 
the mixing angles for three neutrinos almost uniquely PJlOf. Furthermore, analysis in three 
flavor mixing scheme indicates sizeable oscillations of electron neutrinos to tau neutrinos 
that should be observed by the long baseline neutrino experiments such as those utilizing a 
muon storage ring at Fermilab |12[ . These analysis ||[10| also indicate that the solar neutri- 
nos observed on earth should show no MSW effect [flj as the large mass squared differences 
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has been considered in those analysis. Precise measurement of the multi-GeV, 'overhead' 
(cos# 2 ~ 1) events at SuperKamiokande will also be able to verify the three flavor mixing 
scheme |,|TJ as the double ratio R = (%/N e ) measur J(N^/N e ) no osciUation for electron and 



muon for those events is somewhat less than 1 in the three flavor mixing scheme but this 
ratio is 1 in the analysis with a single oscillation process with small mass square differences 
for neutrinos. However present SuperKamiokande data are inconclusive in this low L/E 
region. Three flavor mixing schemes |||10| give very good fit to the SuperKamiokande data 
|| for the double ratio for upward going events ( cos9 z < —0.6 ) but do not give very good 
fit to the data on individual ratio for electron and muon. However, the double ratios are less 
sensitive to systematic errors than the individual ratios. In these analysis f9|J10|l the LSND 
result has been considered as an oscillation effect rather than an unexplained background. 



In near future the BooNE experiment |T3| will test the same channel of neutrino oscillation 
as LSND with higher sensitivity and statistics. Particularly the solutions for the mass square 
differences and the mixing angles in the three flavor mixing scheme as obtained in reference 
H are not significantly contradicted by any existing experimental result and the conflicting 
evidences are below two sigma level. Future various experiments on neutrino oscillation and 
some of those experiments with higher statistics and lesser systematic errors will be able to 
verify the three flavor mixing scheme [pJlOl and it will be certain whether we really need a 



fourth sterile neutrino fL4 |. At present, we feel that three flavor mixing scheme for neutrinos 
are very interesting as it has some specific predictions as mentioned before which can be 
verified by experiments and it tells about the mass squared differences and the mixing angles 
almost uniquely. 



The uniqueness of the mass square differences and the mixing angles [|9].[TI| for three neu 



trinos may have strong impact on physics beyond standard model in the way those constrain 
the parameters of other theories. We like to study such impact on the minimal i?-parity vio- 
lating Supersymmetric model where neutrinos can acquire mass. In supersymmetric models, 
.R-parity was introduced as a matter of convenience to prevent fast proton decay. It is now 
realised that the proton lifetime can be made consistent with experiment without invoking 



discrete i?-parity symmetry. If we do not impose conservation of i?-parity in the model, the 
minimal supersymmetric standard model allows the following B and L violating terms in 
the superpotenial [] 

W = * iik L i L*(E i y + K ]k L l Q 3 (D k ) C + X'r k (w) C (D^ C (D k y (1.1) 

Here L and Q are the lepton and quark doublet superfields, E c is the lepton singlet superfield, 
and U c and D c are the quark singlet superfields and i,j, k are the generation indices. In 
the above, the first two terms are lepton number violating while the third term violates 
baryon number. For the stability of the proton, we assume that only the L - violating 
first two terms in the superpotential is non-zero. One may consider some Z n symmetry to 
remove B -violating term in the superpotential^. As discussed later, L -violating couplings 
give rise to masses for Majorana neutrinos through one loop diagrams as shown in figure 
1, which lead to neutrino oscillation phenomena. In this work we like to show that in 
the R - parity violating Supersymmetric models, it is possible to obtain the required mass 
squared differences and the mixing angles for such massive neutrinos to explain LSND, solar 
and atmospheric neutrino oscillation experiments. In our analysis, it is possible to satisfy 
the bound on the effective mass for the Majorana neutrinos obtained from the neutrinoless 
double beta decay experiment. We have estimated the magnitude of some of the lepton 
number violating couplings and X' i j k which is required to obtain the appropriate mass 
square differences and the mixing angles for neutrino oscillation. This kind of study was 



made earlier [J6| in the two flavor mixing scheme with the lesser available neutrino data. 



1 One may consider another term [i a L a H2 in the superpotential ||31[J . However in general this 
lepton number violating term can be rotated to the first two terms in the superpotential in (1.1) 
unless a symmetry of W does not commute with the SU{4) symmetry of L a rotations in the field 
space. 

2 See [15 1 for other alternative approaches to forbid dimension four as well as dimension five B 



violating operators but keeping L violating operators in the Lagrangian . 



Very recently some other studies 0] also has been made to analyse solar and atmospheric 



neutrino data in the context of Supersymmetric Models. However in our work, unlike other 
works, we have considered solar, atmospheric neutrino oscillation experiments as well as 
LSND data to reconcile with R parity violating Supersymmetric Model. We have also 
discussed the case for which neutrinos may be considered as dark matter candidate. 

There are stringent bounds on different Xijk and \'^ k |1| from low energy processes fl9[ 



and very recently the product of two of such couplings has been constrained significantly 
from the neutrinoless double beta decay |HJ] and from rare leptonic decays of the long-lived 



neutral kaon, the muon and the tau as well as from the mixing of neutral K and B meson 
22fl . In most cases it is found that the upper bound on X' ijk and A^ varies from 1CT 1 to 



10~ 2 for the sfermion mass of order 100 GeV. For higher sfermion masses these values are 
even higher. In our analysis, it seems that for sfermion mass of the order of 100 GeV various 
L violating couplings are less than 10 -2 . Considering the values of L-violating A' couplings 
as obtained from our analysis and also considering the constraint obtained from the electric 
dipole moment of electron it is possible to obtain constraint on the complex phase of some 
Supersymmetry parameters. In section II, we briefly discuss the constraints on masses and 
mixing for three neutrinos obtained from LSND, solar and atmospheric neutrino oscillation 
experiments. In section III, we discuss about masses of three neutrinos in R violating Su- 
persymmetric Model and show how it is possible to reconcile the masses and the mixing 
angles as obtained in the three flavor mixing scheme with the R violating Supersymmetric 
Model. We present the required values of some L violating couplings which satisfy partic- 
ularly various neutrino oscillation experimental data and also satisfy the constraint on the 
effective mass for Majorana neutrinos in the neutrinoless double beta decay experiment. We 
compare these values with the earlier constraint on such couplings. In section IV, we discuss 
that under some circumstances it is possible to get constraint on the complex phase of some 
Supersymmetry parameters like A - parameter. In section V, as concluding remarks we 
mention the possible implications of the obtained values of L violating couplings in collider 
physics and cosmology. 
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II. CONSTRAINT ON NEUTRINO MASSES AND MIXING 



We first mention here the necessary parameters for the three flavor neutrino oscillation. 
After that following references 



we shall consider some specific values for the the masses 
and mixing as solutions to satisfy various available experimental data. The neutrino flavor 
eigenstate are related to the mass eigenstate by 



where U a i are the elements of a unitary mixing matrix U , v a 



e,ii,T 



and Ui 



(2.1) 



^1,2,3- 



According to the standard parametrization [23] of the unitary matrix 



/ 



C12C13 S12C13 Si3^ 

_s 12 c 23 — C12S23S135 C12C23 — Si2S2zSiz5 S23 c 13 
S12S23 - Ci2C2 3 Si 3 (5 -C12S23 - Sl2C23Sl3(5 C23C13 



(2.2) 



where 5 = e l5vi corresponds to the CP-violating phase (which will be neglected here ) and 
c and s stand for sine and cosine of the associated angle placed as subscript. The non 
-diagonal neutrino mass matrix M v in the flavor basis is diagonalised by the unitary matrix 

U v as 



U L V M V U V = D u 



(2.3) 



where D v is the diagonal mass matrix with the real eigenvalues. In the three generation 
neutrino mixing scheme, there are two independent mass square differences. These may be 
considered as A 2 i and A 32 where 



Aij = \m 2 i — m 2 . 



(2.4) 



and rrii and rrij are the neutrino mass eigenvalues. From the solar neutrino deficit one may 
consider M 



10~ 4 eV 2 < A 21 < l(T 3 eV 



(2.5) 



in which the lower limit is obtained from SuperKamiokande data and the upper limit 
is obtained from the CHOOZ experiment |nj| . Keeping in mind both the atmospheric and 
LSND data, another mass square difference A32 can be considered as || 

A 32 « 0.3eV 2 (2.6) 



in which the lower limit from the Bugey reactor constraint |2J] and the upper limit from the 



CDHSW |25[ have also been considered. In the one mass square difference dominating the 
other, the three flavor mixing scheme greatly simplifies and one can write the probability of 
the observation of neutrino oscillation in LSND and SuperKamiokande in terms of A32 and 
the elements of U . For the solar neutrino experiment the sine squared terms containing 
two mass square differences in the expression for the survival probability of solar electron 
neutrino, can be averaged for the flight length and the energy of the neutrinos observed on 
earth. In this case the probability can be written in terms of the elements of U only. As 
the probability of oscillation in LSND and Superkamiokande and the survival probability for 
solar electron neutrinos are provided by the experiments, one can solve for the three angles 
by which matrix U in (2.2) is defined. The results obtained by Barenboim et al |J show 
that four set of solutions for three angles are possible. However, two set of solutions can 
be discarded by considering the SuperKamikande zenith angle (cos #.2 < —0.6) behavior of 
atmospheric neutrino data for upward going events. The other two allowed set of solutions 
for the three angles as obtained in reference |J are 

12 = 54.5; ^13 = 13.1; 6 23 = 27.3 (2.7) 



#12 = 35.5; ^13 = 13.1; 6 23 = 27.3 (2.8) 

The result obtained by Thun et al [[1(J to satisfy solar , atmospheric and LSND data matches 
almost with the second set of solutions for three angles as mentioned in (2.8). In our analysis 
we shall consider either (2.7) or (2.8) for the three angles which specify the unitary matrix 
U in (2.2). 
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The neutrino oscillation experiments give us information about the mass squared differ- 
ences of three neutrinos in the three flavor mixing scheme as discussed above. However, to 
know the mass of different neutrinos we have to consider some other experiment. The masses 
are generated for Majorana neutrinos in R violating Minimal Supersymmetric Model, so we 
have to consider the constraint coming from neutrinoless double beta decay. This gives us 
an estimate for the masses of neutrinos. The contribution of Majorana neutrinos to the 
amplitude of the neutrinoless double beta decay j32 



is 



| M | = | J2 U ei 2 m Ul |< m{0v(3/3) = 0.68 eV 

i=l,2,3 

To satisfy this constraint and keeping in mind that there are some uncertainties in the 
calculation of the nuclear matrix elements one may consider different masses of the Majorana 



neutrinos of the order of eV or less p9 |. Considering (2.5) and (2.6) alongwith this constraint 
it is found that there are two interesting possibilities for the masses of neutrinos. In one 
case, all three neutrinos have almost degenerate mass and we may consider 

m 2 « leV (2.9) 

then the masses for other two neutrinos are 

mi « leV; m 3 « l.UeV (2.10) 

In another case, the masses of two neutrinos are nearly degenerate whereas the third one is 
heavier and we may consider 

m 2 ^3xlO" 2 eV (2.11) 

then the masses for other two neutrinos are 

mi » 2 x 10~ 2 eV; m 3 » 0.55eV (2.12) 

One may consider the neutrinos as candidate for the dark matter solutions also. In that 
case, if one assumes Q = 1 and the energy density of the neutrinos p u = 0.2p c where p c is the 
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critical density in the Big-Bang Model p8| , it is desirable to have the sum of the neutrino 
masses around 5 eV and one may consider the nearly degenerate three masses of neutrinos 
given by (2.9) and (2.10). 

In our analysis we shall consider the above mentioned four interesting possible solutions 
for masses and mixing angles - one set of solutions from (2.7), (2.9) and (2.10), one set of 
solutions from (2.8), (2.9) and (2.10), one set from (2.7), (2.11) and (2.12) and the other set 
from (2.8), (2.11) and (2.12). We shall discuss in the next section how all these solutions 
can be reconciled with R parity violating Supersymmetric Model. 

III. NEUTRINO MASS MATRIX IN R- VIOLATING SUPERSYMMETRIC 
MODEL AND CONSTRAINT ON L VIOLATING COUPLINGS 

The trilinear lepton number violating renormalizable term in the superpotential in (1.1) 
generates Majorana neutrino masses |]30| , |3T| through the generic one loop diagram as shown 



in figure 1 in which s and s stand for either lepton and slepton or quark and squark respec- 
tively. The helicity flip on the internal fermion line is necessary and that requires the mixing 
of s and s c . The contribution to the mass insertion as shown in figure 1, is proportional 
to the mass of the fermionic superpartner of Sk and is also proportional to m ( ~ A 
~ /i, where A and \i are the susy- breaking mass parameter) . The single diagram in figure 
1 that contributes to the Majorana neutrino mass matrix m UiVj is 

m ViV ss - 2 3.1 

when one considers the lepton and slepton for s and s in the diagram in figure 1. Both the 
diagrams in (a) and (b) are to be considered together and summed to evaluate the neutrino 
mass matrix element. However for i — j and k = n, the two diagrams coincide and for that 
only one is to be considered. For quark and squark in the diagram the similar contribution 
will be obtained. However, in that case, the above contribution is to be multiplied by a 
color factor 3 and fhk in the above equation is to be considered as the squark mass instead 
of slepton mass and the A couplings in (3.1) is to be replaced by A' couplings. 
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In constructing neutrino mass matrix we shall consider the following things. Firstly we 
shall relate squark and slepton mass as 

^slepton = ^squark/ K (3-2) 

where K is a number depending on the various choices of Supersymmetry parameters. Differ- 
ent squarks have almost degenerate mass and different sleptons also have almost degenerate 
mass as otherwise there is severe constraint from the flavor changing neutral current. There 
are 9 A couplings and 27 A' couplings entering the neutrino mass matrix. However, in writing 
each element of the neutrino mass matrix we shall consider only the leading term in terms 
of the magnitude of mass obtained from (3.1). We shall consider the diagram with lepton 
and slepton in figure 1 and shall also consider the diagram with squark and quark in figure 
1 in each element of the neutrino mass matrix for which two different types of L violating 
couplings appear in each element. Under this consideration only the following L violating 
couplings appear in the neutrino mass matrix. These are A' 133 , A' 233 , A 333 , Ai 33 , A233, A232, Ai 32 
couplings. The notations for the first five couplings in later discussion will be A^ , A^, A^ , 
Ag , A^ respectively. We are ignoring the effect of other couplings in our analysis and we are 
assuming that the A couplings are not much hierarchical among themselves and A' couplings 
are also not much hierarchical among themselves. At the end of this section we shall make a 
few qualitative comments about considering other L violating couplings in the mass matrix. 
We write the neutrino mass matrix as 



N = a 



Km 2 T X l e 2 + 3m 2 \f 2Km 2 X l e X l lx + 6ml XIX^ -IKm^m^X^ + 6mgA«A«^ 

2Kml\% + 6m 2 A«A« Km\\^ + 1m\\f -2Km T m^X 232 + 6m 2 A« A* 



-2Km^m T X l ^X 132 + 6mgAfA« -2Km T m^\ 232 + 6m fe 2 A£ A« Km\ \ 2 232 + 3m 2 b \f 



J 

(3.3) 



where 



Tfl 

a = i6^l (3 - 4) 
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and m s is the almost degenerate squark mass. The eigenvalues for this matrix correspond 
to three masses m 1; m 2 and m 3 for three Majorana neutrinos. We can write the diagonal 
mass matrix as 



^mi 



D„ 



(3.5) 



/ 



m 2 
^ m 3 

All the elements of this diagonal mass matrix can be written by considering particular set of 
solutions for the masses from the earlier section. The unitary matrix U u in (2.2) diagonalising 
the non-diagonal neutrino mass matrix in the flavor basis is also known to us if we consider 
particular set of solutions for the three angles from the earlier section. As both D u and U v 
are known we can obtain the non-diagonal mass matrix M v in the flavor basis using the 
relation 



U v D v Ul = M v . 



(3.6) 



So for particular set of solutions for the masses and the mixing angles discussed in the earlier 
section, all the elements of M v are known. However, this M v is equal to N which is also the 
non-diagonal mass matrix expressed in terms of different L violating couplings. So we write 



N — M,. 



From (3.7) we get six equations for the L violating couplings 



Km 2 T X l e 2 + 3m 2 b \f = M„(l, I) /a 



(3.7) 



(3.8) 



2Km 2 T X l e X l tl + 6m^A^ = M v (l,2)/t 



2Km IJ ,m T X l lx X 13 2 + Qm 2 b X q e X q T = M„(l, 3)/a 



Km 2 X 1 ^ + 3ml\f = M„(2, 2) /a 



(3.9) 
(3.10) 
(3.11) 
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- 2^771^^232 + 6m*A«A« = M u (2, 3)/o 



(3.12) 



Km l^ 2 232 + 3mgA« 2 = M„(3, 3)/a (3.13) 

after comparing the elements (1,1), (1,2), (1,3), (2,2), (2,3) and (3,3) respectively. However, 
there are seven L violating couplings involved in these six equations. So we shall consider 
the possible value of one of the L violating couplings from some other experiment instead 
of neutrino oscillation experiments for solving above six equations to find six L violating 
couplings. We shall consider particularly some value of A232 lower than °'°^ ls which is 
allowed after considering the constraint from lepton universality fT8||T9|j . From these six 



equations we can determine the values of six L violating couplings for which it is possible 
to reconcile LSND, solar and atomspheric neutrino oscillation experimental data with the 
R parity violating Supersymmetric Model. 

To determine M v we first consider (2.7), (2.9) and (2.10) for the masses and the mixing 
angles. From (2.9) and (2.10) we get a specific D v in (3.5), and from (2.7) we get a specific 
U v in (2.2). Using relation (3.6) we obtain the following form of M v ( Here and in later 
discussions to obtain M v we shall consider m 1; m 2 and in (3.5) in eV unit) : 



My 



( 1.00712 0.0143034 0.0274612^ 



0.0143034 1.02782 0.0542459 



1 0.0274612 0.0542459 1.10499 j 



(3.14) 



We take m b = 4.3 x 10 9 eV , m T = 1.777 x 10 9 eV and m M = 0.105658 x 10 9 eV and 
solve (3.8)-(3.13) after considering a specific M v in (3.14). For various allowed real values 
of A232 lower than that mentioned earlier, the solution for other six L violating couplings 
do not change by an order. We present below the the values of these couplings considering 
A232 in the range (10~ 6 — 10~ 3 ) ^= (Here and in later discussions m s and fh stand for the 
corresponding magnitude in GeV unit): 

5.3xl0- 5 m s 1.3 x 10^m s 5.6 x 10~ 5 m s 
K « 7^ > K ~ 7^ > K ~ 7^ ( 3 - 15 ) 
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, 7.3xl0- 6 m s , 2.3xl0- 4 m s 4.0 x 10- 3 m s 

e ~ 7^ — ' — ' 132 ~ 7^ — ^ ' 

V Km v Km v Km 
Another set of real solutions for various L violating couplings for the above case is given 
below : 

5.3xl0- 5 m s 1.3xl0- 6 m s 5.6 x l(T 5 m s 

^ « 7 = , A« « 7 = , A* « 7 = (3.17) 



, 4.2xl0- 6 m s , 2.3xl0~ 4 m s 3.9 x l(T 3 m s 

A e ~ , , A ~ . , A132 ~ , ^o.loj 

yKm yKm ' y Km 

We have ignored overall + or - sign for the solutions (here and in later cases also) for different 
L violating couplings. Although there are different set of solutions possible but if we ignore 
the small changes in the higher decimal places for different solutions then mainly the two 
set of solutions are found to differ to some extent from each other particularly for the value 
of A e and A132 and those two sets of solutions are presented above. 

It is important to note here that there is almost no change of the A132 value for various 
real A232 value in the range (0.0 — 10~ 3 ) -7= and the two values of Ai 32 in (3.16) and (3.18) 
are very near to the upper bound obtained from the experimental value of R T = f?^zf~^ 
|T8^PU| . This indicates that there is possibility to see the lepton universality violation in fu- 
ture experiments. Same comment is also true for A232 coupling as almost same real solutions 
for various L violating couplings exist for the higher allowed value of A232 coupling also. 
Of course the statement is based on the present neutrino oscillation data and considering 
neutrino as Majorana particle , the main contribution in the neutrino mass matrix coming 
from the earlier mentioned seven L violating couplings and the L violating couplings consid- 
ered here being real. Although in obtaining the solutions for L violating couplings we have 
considered here almost degenerate mass for three neutrinos , but such higher values of A 132 
or A232 are possible for hierarchical nature of the masses of neutrinos also, as can be seen in 
the later part of our analysis. If we consider complex or imaginary value of A232 coupling 
considering the experimental upper bound mentioned earlier it is possible to obtain complex 
solutions for other six L violating couplings from (3.8)-(3.13) for M v in (3.14). For brevity, 
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we are not presenting those solutions of various L violating couplings for this case of masses 
and mixing angles and for other cases also. However, at the end of this section we shall make 
a few general remarks on the complex solutions for these L violating couplings. Existence 
of the possible solutions for L violating couplings in (3.15) and (3.16) or (3.17) and (3.18) 
indicates that considering the almost degenerate mass neutrinos (which may be candidate 
for dark matter also) as mentioned in (2.9) and (2.10), it is possible to reconcile LSND, 
solar and atmospheric neutrino oscillation data with the R parity violating Supersymmetric 
Model. 

Next, we consider the other possible solutions for the mixing angles as stated in (2.8) and 
consider again the almost degenerate mass of neutrinos as mentioned in (2.9) and (2.10). In 
this case, we get the following form of M v after using (3.6) : 



1.00704 0.0143116 0.0274772^ 



M„ 



0.0143116 1.02788 0.054211 
0.0274772 0.054211 1.105 



(3.19) 



Like earlier case, we again solve equations (3.8)-(3.13) for specific M v in (3.19) and obtain 
the solutions for six L violating couplings. All the solutions in this case are approximately 
same as (3.15)-(3.18). So the different set of choices for the mixing angles in (2.8) do not 
lead to significant change in the values of L violating couplings. 

We shall consider next the hierarchical neutrino masses as mentioned in (2.11) and (2.12). 
For the three mixing angles we consider (2.7). As before using (3.6) we obtain the following 
form of M v : 



^0.0534391 0.0569347 0.10027 X 



0.0569347 0.127333 0.202493 
0.10027 0.202493 0.417771 



(3.20) 



Solving (3.8)-(3.13) for specific M u in (3.20) we obtain the following real solutions for six L 
violating couplings. We present below the values of these couplings considering A232 in the 
range (10~ 6 - 10~ 3 ) jk. 
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Another set of real solutions for various L violating couplings for the above case is given 
below : 

i.3 x 10- 6 m, .„ 3.4xl0~ 5 m, 



1.2 x l(T 5 m s 

i - 



V 



A? 



(3.23) 



A' « 



5.2 x 10" 6 m s 



A 



7.2 x 10" 5 m s 



A 



1.2 x 10" 3 m s 



\jKvn 11 \jKm 132 V 'Km 

So it is seen that considering hierarchical nature of the masses of neutrinos also it is possi- 
ble to reconcile LSND, solar and atmospheric neutrino oscillation experiments with the R 
violating Supersymmetric Model. In this case only thing to note here is that A132 is slightly 
lower than the earlier cases however not far from the present experimental bound obtained 
from the lepton universality violation [fX8| |2Q | . 

Next we shall consider the hierarchical mass pattern of neutrinos like earlier case but 
consider the mixing angles as presented in (2.8). In that case, using (3.6) we get the following 
form of M„ : 



(3.24) 



^0.0503506 0.0572644 0.100908^ 



M,, 



0.0572644 0.129869 0.201098 
0.100908 0.201098 0.418324 



(3.25) 



Like before we solve (3.8)-(3.13) for M u in (3.25) and consider the same range for A232 like 
earlier cases. In this case, the solutions for L violating couplings are almost same as before 
with hierarchical masses of neutrinos and we are not presenting those solutions seperately. 

If the future neutrino oscillation experiments with higher sensitivity and more data 
support the three flavor mixing scheme as mentioned in section II and the L violating 
couplings are real, it is expected that experiments on lepton universality violation in future 
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will find signal for the values of A132 or A232 couplings at the level required by our analysis. 
However, if no signals are found for those values of L violating ocuplings- particularly A132 
coupling then the explaination for that may be the following. In that case normally it will 
be expected that A132 coupling is very small. However then if one considers again those 
six equations (3.8)-(3.13) considering Ai 32 as effectively zero it is found that for the various 
cases for masses and mixing angles, real values of A232 has to be always of the order of 
W yjKih ' However if the signal for A232 is also not seen through r- universality violation, it 
will be necessary to check the role of other couplings for the analysis of neutrino masses and 
mixing angles. As under this circumstance, A132 and A232 will be smaller, we may consider 
the terms next to the leading order in mass in the various elements of the neutrino mass 
matrix in (3.3). As the mass factor associated with those other non-leading contribution 
will be less in magnitude, it is expected that the magnitude of some other coupling should 
be somewhat higher like A132 to reproduce the similar forms of M v mentioned earlier and the 
lepton number violation, in that case, should be observed through that coupling. Depending 
on the results of future experiments on neutrino oscillation, tau universality violation etc., 
the analysis with other such couplings may be important. 

We like to make a few remarks on the complex solutions for various L violating couplings. 
Earlier in obtaining all the above-mentioned solutions for different L - violating couplings 
we have considered value of A232 coupling in the range (10~ 6 — 10~ 3 ) J^^ - However, if 
one considers the value of A232 in the range (4.0 — 6.0) 10~ 3 which is very near to the 



experimental upper bound fl8 |, then from the equations (3.8)-(3.13) considering different 
form of M v as mentioned earlier, one will obtain the complex solutions for various L violating 
couplings. Furthermore, if one considers imaginary or complex value of A232 in that case also 
one will obtain the complex solutions for various L violating couplings. Depending on the 
various choices of the values of A232 one may obtain from (3.8)-(3.16) the various solutions 
for different L violating couplings with various possible complex phases. For brevity, we have 
not presented various possible complex solutions. For some complex phases the solutions 
may not be allowed depending on constraint particularly from the value of the electric dipole 
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moment of electron. We have discussed this in the next section. 

The solutions for all A and A' couplings are obtained in terms of the parameters rh, 
rh s and K. Here rh is the Supersymmetry breaking mass parameter which is expected 
to lie in the range of 0(100 GeV ) to 0(TeV ). To compare our constraint on A and A' 
couplings with the earlier constraints JT^] we shall consider K pa 1 which means slepton 
mass does not differ much from the squark mass and also shall consider both squark mass 
and rh of 0(100 GeV ). We shall consider those solutions of L violating couplings for which 
complex phases are negligible. The earlier constraint from the tau universality violation 
is Ai 32 < 0.06, A 232 < 0.06 and A 233 = Aj, < 0.06 which in our analysis from neutino 
oscillation, are found to be 0(0.01-0.04), 0(0.01) Q and 0( 0.0007-0.0023) respectively. The 
earlier constraint from R t = T hadron {Z°) /Ti{Z°) @ A' 333 = 1% < 0.26 and A' 233 = l%< 0.39 
which in our analysis are 0(0.0003 — 0.0005) and 0((1 — 8) x 10~ 5 ) respectively. The earlier 
constraint on A' 133 = \ q e obtained from the constraint from neutrino mass ]34|] is \ q e ~ 0.002 
and A 133 = X l e pa 0.004 which in our analysis are 0(0.00008-0.00053) and 0(0.00004-0.00037) 
respectively. So our analysis indicates somewhat lower values of various L violating couplings 
than the upper bound on these couplings obtained from other experiments. Furthermore, if 
one considers rh to be nearer to TeV region then these values of L violating couplings will 
be further lowered. The upper bounds on these couplings obtained from the experimental 
data on neutral currents , (3 decay fL9"| , muon decay (fi — > ej,fi — > eee) fl35"| , or tau decay 
(r — > /i7, r — > ej) |36[] etc. are somewhat higher than the values required in our analysis. 



3 When we consider very small value of A132 which can be neglected in the neutrino mass matrix in 
(3.3), in that case we get this solution for A232 from the equations (3.8)-(3.13). Otherwise various 
solutions for X232 are possible as mentioned earlier. 
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IV. CONSTRAINT ON COMPLEX PHASE OF SUPERSYMMETRY 

BREAKING PARAMETER 



In the Standard Model the electric dipole moment of electron is much smaller than their 
present experimental bound d e < 10 _26 ecm [|37|. So the new sources of CP violation which 



occurs in the supersymmetric model can be studied on the basis of electric dipole moment 
of electron P%|-4*U|. In the minimal supersymmetric standard model apart from the Yukawa 



couplings there are several complex parameters like three gaugino masses corresponding to 
SU(3), SU(2) and U(l) groups, the mass parameter m# in the bilinear term in the Higgs 
superfields in the superpotential , dimensionless parameters A and B in the trilinear and 
the bilinear terms of the scalar fields. With suitable redifinition of the fields, some of 
these parameters can be made real but in that case some others can not be made real like 
A parameter |39fl . The complex A will contribute to the electric dipole moment (edm) of 
electron. Furthermore, if we consider the complex A' couplings, the complex phase associated 
with those will also contribute to edm of electron. There will be various diagrams in the im- 
parity violating Supersymmetry for the edm of electron |40| . But the significant contribution 



to edm comes from the one loop diagram containing top quark in the loop as shown in Figure- 
2. There will be diagram containing massive neutrinos in the loop . However, the masses 
of neutrinos are quite small in our discussion and we are ignoring those types of diagrams 
for our discussion as there will be lesser contribution to the edm of electron. In terms of 
complex phases we can write Af and \'^ k as 

A u4 =| A | exp (ia Aud ), X' ijk =| X' ijk | exp (4.1) 

and the mixing angle for the left and the right squark in the familiar way as 

tan 26 = 2 | A u \ m u /(/i 2 L - fi 2 R ). (4.2) 



Following reference j4(J and assuming different \[ fc containing another complex phase as 



mentioned in (4.1) we can write the edm of electron from figure 2 as 
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d e « — sin^ cos#(^(cos /3 — sin (3) sinaA d + cos/3 sin/3 cos«a 

x | A' life | 2 2 ^ 1 m Mj [Fi(a*) + 2F 2 (x fc )] 10- 17 e cm (4.3) 

where Xk = (md k /m s ) 2 and the loop integrals F\ and F 2 are expressed in terms of Xk as 

n / \ 1 / 2xfc lnxfc\ 

2(1 -x fc ) V / 

*!.(**) = - 1 2 (s-x k + . (4.4) 

2(1 -x k ) V 1-^fc/ 

Ad and m s in (4.3) correspond to the magnitude of those quantities expressed in GeV. We 
are particularly interested for j = 3 and k = 3 case in (4.3). We got a solution for \ q e = A' 133 
in section III to explain the neutrino physics data. So we like to constrain here particularly 
the complex phases associated with A in (4.3). Considering as real or complex and writing 
it as Cm s j Vm in the form obtained in section III (where C is some value depending on the 
type of solutions), A and m s both from 100 GeV to 1 TeV it is found that the constraint on 
the complex phase ctA d and (3 is 



cos 2 (3 — sin 2 0) sin«A d +cos/3 sin/3 cosaA d J sin 6* cos# 

< (2-16)xl0-» (4 5) 

In the case for which there will be no contribution to edm in (4.3) is f3 = —a a/2- For 
the complex solutions of \ q e for which | C \ is not less than about 10~ 6 we can make the 
following statements. For (3 ~ 7r/4 and a.A d ~ 7r/2 the above inequality can be satisfied for 
any value of 9. However such a large phase for A is not possible as the edm of electron will 
get contribution from other diagrams involving charginos and neutralinos at the one loop 
level p8| cancelling this possibility. So (3 ~ 7r/4 is not possible for any value of 9. So those 



set of complex solutions for different L violating couplings should not be considered when 
the complex phase associated with \ e q is found to be approximately 7r/4 and | C \ satisfies 
the above condition. Let us consider that A' ljfc are real and 9 = 7r/4 in that case it is seen 
from section III that C ~ 10 -5 and the complex phase for the A parameter otA d < 3.2 x 10~ 3 . 
Without any specific choice of the mixing angle 9 one can constrain only the combination 
of (3, a a and 9 as shown in (4.5). 
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V. CONCLUSION 



We have shown here that in the minimal supersymmetric model with .R-parity violating 
trilinear term in the superpotential in (1.1) it is possible to obtain the appropriate mass 
square differences and the mixing angles as required to explain the LSND, atmospheric and 
solar neutrino oscillation experimental data in the three flavor mixing scheme for neutrinos. 
The validity of the three flavor mixing scheme can be verified in the near future experiments 
on neutrino oscillation as mentioned in the introduction. The masses for three Majorana 
neutrinos are generated at the one loop level as shown in Figure-1 and it is possible to satisfy 
the constraint on the masses and the mixing angles coming from the neutrinoless double beta 
decay. In each element of the neutrino mass matrix in (3.3) we have considered two leading 
terms in terms of the magnitude of masses in (3.1) coming from the diagram with slepton 
and lepton and also coming from the diagram with quark and squark in Figure-1. Under 
this consideration it is interesting to note that for real values of various L violating couplings 
at least one of the couplings either Ai 32 or A232 is expected to be quite high and very near 
to the experimental upper bound coming from the r -universality violation [p~8| pOf| ■ Apart 



from these two particular couplings for some of the L violating couplings the magnitude are 
such that it might be possible to observe such L violating interaction at the Tevatron or 
at HERA. At the Tevatron after squark pair production those squarks will decay to LSP 



(say neutralino) and which will decay via UDE operator giving multilepton signal [41 



At HERA one can see R violating Supersymmetry signal for L % Q^D kC operator [42] through 



resonant squark production and its subsequent decay to electron or positron and neutrino 
giving the signal of high p? electron or high positron or missing px for neutrino. The 
basic requirement for the observation of such signal is that LSP has to decay inside the 
detector and this puts bound f4l||T2| 



UOOGeW 

where mj- stand for the squark and the slepton mass. From the above condition it is 
seen that if we consider fh of the order of squark mass fh s in that case for squark mass or 
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slepton mass of the order of 300 GeV, it may be possible to observe L violating signal for 
those couplings discussed in section III for which A, A' > 5 x 10 -6 -^; whereas for squark 
or slepton mass of the order of TeV the condition is A, A' > 3 x 10~ 5 ^=. So for various 
couplings considered in our analysis in section III, it might be possible to observe L violating 
signal. 

If one considers the baryogenesis in the early universe at the GUT scale, after the gen- 
eration of asymmetry to satisfy the out of equilibrium condition one requires L violating 
couplings ~ 10~ 7 for squark mass from 100 GeV to 1 TeV range [43| which is significantly 
smaller than the values of some of the couplings obtained in section III. So if one likes to 
satisfy the neutrino physics experimental data in the three flavor mixing scheme, it seems 
in the R violating Supersymmetric scenario the generation of the baryonic asymmetry near 
the electroweak scale is more favored where the constraint on L violating couplings are not 



so severe 



44]. We have shown that in the R violating Supersymmetric models neutrino 



can be considered as dark matter candidate also. Our analysis also indicates that to satisfy 
various experimental data on neutrino oscillation , the lepton number violating couplings are 
constrained in such a way that some combinations of left and right squark mixing angles and 
the complex phases of some Supersymmetry parameters - particularly that of A parameter 
are constrained. 
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FIG. 1. One loop diagram involving L-violating couplings generating neutrino mass. 
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FIG. 2. One loop diagram contributing to electric dipole moment of electron. 
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